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Abstract. We discuss the structure of positive definite kernels in terms of 
operator models. In particular, we introduce two models, one of Hessenberg 
type and another one that we call near tridiagonal. These models produce 
parametrizations of the kernels and we describe the combinatorial nature of 
these parametrizations in terms of lattice paths of Dyck and Lukasiewicz type. 



1. Introduction 

In this paper we are interested in parametrizations and combinatorial descrip- 
tions of positive definite kernels on the set No of non-negative integers. Positive 
definite kernels are complex valued maps K : No X No — > C with the property that 
for each n > and each choice of elements pi, . . ., p n in No and complex numbers 
Ai, . . ., A„, we have 

n 

(1.1) J2 K{p k , Pj )\^\ k >Q. 

A fundamental result of Kolmogorov, |S] provides a Hilbert space interpretation of 
positive definite kernels as Gram kernels, that is, there exists a Hilbert space 7i 
and elements v(n) 6 H, n > 0, such that 

(1.2) K(i,j) = (v(j),v(i}). 

Two of the best known examples of positive definite kernels are those of Toeplitz 
type, for which K(i + l,j + I) = K(i,j), I S No, and those of Hankel type, for 
which K(i,j + I) = K(i + l,j), S No- In both these cases the representation 
(|1.2|) can be improved by some more specialized descriptions that might be called 
operator models of the kernels. Thus, if K is Toeplitz (for simplicity we assume 
K(0, 0) = 1 and all the inequalities in 1|1.1[) are strict) then there exists an isometric 
operator W, written in upper Hessenberg form, such that 



(1-3) K(i,j) = e W j - 



: o- 



where eo = [ 1 . . . ] . Likewise, if K is Hankel (with same simplifications as 
above) then there exists a symmetric operator J, written in tridiagonal form, such 
that 

(1.4) K(i,j) = e J i+ ieZ, 

Our goal is to extend both models i|1.3fl and i|1.4fl to arbitrary positive definite 
kernels on No without Toeplitz or Hankel assumptions. These models will produce 
parametrizations of the kernels and we will give combinatorial descriptions of these 
parametrizations in terms of lattice paths. 

l 



2. Isometric Hessenberg models and Dyck paths 



In this section we show that any positive definite kernel on No has a Hessenberg 
model and then we show how to relate this model to the set of Dyck paths. In 
order to simplify the notation we consider only positive definite kernels for which 
all the inequalities in Ijl.ljl are strict, when we say that the kernel is stricly positive 
definite, and we also assume K(l, I) = 1 for all I > 0. Both these assumptions 
can be easily removed. In addition, all our considerations can be easily adapted 
to kernels K : No X No — > £(£), where £(£) denotes the set of linear bounded 
operators on the Hilbert space £. 

We now introduce the elements necesary in the presentation of the results. For 
a complex number 7 with I7I < 1 we define its defect by d 1 = (1 — I7I 2 ) 1 / 2 and its 
Julia matrix by 

7 d 7 
d~t -7 

We note that the Julia matrix is unitary and this construction can be extended to 
certain families of complex numbers as follows. Let Y = {"fk.j}o<k<j be a family 
of complex numbers such that |7fe,j| < 1 for all k < j. For simplicity we will write 
dkj instead of d lh j . We can now describe the Hessenberg model. First we define 
for k < j, 

(il © J{lk,k+2) © Ij-k-2) ■ ■ ■ (Ij-k-1 © J(lk,j)) ) 

where /; denotes the I x I identity matrix. Then we introduce the operators Wk(T) 
on the Hilbert space Z 2 (No) of square-summable sequences by the formula: 

W k (T) = s - lim (V k ,j(T) e 0), k > 0, 

j->oo 

where s — lim denotes the strong operator limit. It is easily seen that each Wfc(r) 
is an isometry with upper Hessenberg matrix with respect to the standard basis of 
Z 2 (N ), that is, if (W fe (r)) iJ denotes the entry of W k (T), then (W k (T))ij = 
for i > j + 1. It is also useful to consider the unitary matrices Uk,j(T) defined 
recursively by J7fc,fc(r) = 1\ and for k < j, 

U k , j (r) = V k , j (T)(U k+hj (T)®I 1 ). 

We can prove now the existence of an isometric Hessenberg model for any strictly 
positive definite kernel on No- 

Theorem 2.1. If K is a strictly positive definite kernel on No with K(l,l) = 1 
for all I > 0, then there exists a family {W k }k>o of isometric Hessenberg operators 
such that for j > i, 

(2.1) K(i,j) = e Wi(T)W i+1 (r) . . . Wj_i(T)e5. 

Proof. This is just a reformulation of Theorem 2.3 in (see also 0, Chapter 
1). Thus, by Theorem 1.3 in |2j, there exists a uniquely determined family T — 
{lk,j}o<k<j of complex numbers such that K(i,j) — (Uij(T)) Q for i < j. Then it 
is easily seen from the definitions that 

e W i (T)W i+1 (T)...W j - 1 (T)e* = {U id (T)) 0t0 . 

□ 



When if is a Toeplitz kernel, then (|2.1|) reduces to l|1.3|l and the parameters 
7fcj satisfy = 7i,j for i < j, I > 1. The numbers j n = j ,n, n > 1, are 

called the Szego parameters of if (other names, like Schur parameters, reflection 
coefficients, or Verblunski parameters are currently used in the literature), and they 
play a central role in the theory of orthogonal polynomials on the unit circle and its 
many applications, [H] and, for a recent account, [HJ (which also contains a detailed 
discussion of the Hessenberg model in the Toeplitz case). 

Next we explain the connection between Hessenberg models and Dyck paths. A 
Dyck path of length 2k is a path in the positive quadrant of the lattice Z 2 which 
starts at (0,0), ends at (2k, 0), and consists of rise steps / and fall steps \ (see 
Figure 1). For more information on Dyck paths and their combinatorics, see 0. 




FIGURE 1 . A Dyck path of length 8 and height 3 

Let T>k be the set of Dyck paths of length 2k and let Ah be the set of points 
(l,q), q > 0, with the property that there exists p <E T>u with (l,q) € p. It is seen 
that 

A k = {{j + i,j-i) \0<i<j<k}. 
Also, we notice that if p S T>k and x = (l,q) G p, then there are only four types 
of behaviour of p about x: (I) a rise step followed by a fall step; (II) a fall step 
followed by a rise step; (III) two consecutive rise steps; (IV) two consecutive fall 
steps (see Figure 2). 
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Figure 2 . Behaviour of a Dyck path about a vertex x 6 Ak 
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Let p be a Dyck path in Vk with the property that (21, 0) £ p. The restriction of p 
from (21,0) to (2fc,0) is called a Dyck subpath starting at (21,0) in T>k- The set of 
all possible Dyck subpaths starting at (21, 0) in T>k is denoted by V l k and there exists 
a bijection between T> k and T>k-i- 



This implies that the number of elements in T> l k 



is given by the Catalan number Ck-i 



1 



k-l 

3 



1 



2(k - I) 
k - I 



also. 



V k . If 



q G Pi then there could be many Dyck paths whose restrictions at (21,0) coincide 
with q, however if pi and P2 are two such Dyck paths then ajj(pi) = ajj(p2) for 
j + i > 21. We will write ajj(q) in order to denote this common value. 

We now describe the structure of the strictly positive definite kernels on No- 

Theorem 2.2. The kernel K on No with K(l,l) — 1 for all I is strictly positive 
definite if and only if there is a family {"fk,j}o<k<j of complex numbers, \"fk,j\ < 1 
for all k < j , such that 

(2.2) K(l,m)= II a M"(q)- 

qgD^ l<i<j<m 

Proof. Half of this result was proved in [Q, but we give some details here for com- 
pleteness. Assume that K is strictly positive definite. By Theorem 1.3 and Theo- 
rem 2.3 in [5] there exists a uniquely determined family T — {jk,j}o<k<j of complex 
numbers such that K(l,m) = (E/i, m (r))o,Oj the (0,0) entry of the matrix U^ m (T). 
It is convenient to visualize this relation by means of a so-called transmission line, 
as showed in Figure 3 for K(0, 2) and K(0, 3). 




FIGURE 3. Transmission line for K(Q, 3) 



Thus, if 1 is the input at A then at B we read off the expression of K(0, 3) in terms 
of the parameters 7 0) i, 70,2, 7o,3> 7i,2> 7i,3> 72,3 and their defects, 

K(0, 3) = 70,171,272,3 + 70,1^1,271,3^2,3 + 1^0,170,2^1,272,3 
— ^0,l70,27l,27l,3^2,3 + ^0,1^0, 270,3^1, 3^2,3- 

Likewise, if the input at C is 1, then the output at B is now the expression of 
K(0,2), 

K(0, 2) = 70,171,2 + ^0,170,2^1,2 

(for more details see |3]). Each path in the transmission line contributes an additive 
term in K(l, m). Going from a path in the transmission line to a Dyck path is easy, 
each box associated with a Julia matrix corresponds to a point in Ak , see Figure 4. 
It is also clear that each additive term in K(l, m) is given by Y\.i<i<j<m a «,i( f l) f° r 
some q £ V l m . This gives $T$ . 

Conversely, given a family T = {^k,j}o<k<j of complex numbers with |7fc,j| < 1 
for all k < j, we define 

K{l,m) = II a M(q)- 

qeT>' m l<i<j<m 
4 



By the first part of the proof, this gives K{1, m) = (U^ m (T))o^, and it remains to 
show that K is a strictly positive definite kernel on No- By Theorem 2.1, 

K{i,j) = e Wi(T)W i+1 (T) . . . ^_i(r)eS 

for i < j. This relation implies that if is a positive definite kernel. Also, by 
Proposition 1.7 in 

det[if(Z,m)]£ ro=0 = [] d h >0 > 

0<i<j<n 

so that if is a strictly positive definite kernel on Nq. □ 
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Figure 4. From a path in a transmission line to a Dyck path 



Remarks (a) It is quite simple to remove the two restrictions on K considered 
in Theorem 12.21 First, formula (|2.2I) still provides a one-to-one correspondence 
between the set of positive definite kernels on No with K(l, = 1 for all I and the 
set <S of families {"fk,j}o<k<j of complex numbers with the properties: \^k.j\ < 1 
for all k < j; if \^k,j\ = 1 for some pair (k,j), then 7^ = for I < k and 7fc iTO = 
for m > j. 

(b) If we remove the assumption that K(l, I) = 1 for all I, then the diagonal elements 
K (I, I) of the kernel K could be considered as parameters and there is a one-to-one 
correspondence between the positive definite kernels on No and the set 5+ of pairs 
({/;}i>o, {7fej}o<fc<i), where /; > for alU > and {^k,j}o<k<j is an element of 
S with the additional property that if // = for some I > 0, then jk,i — and 
li.m = for k < I and m > I. Formula (|2.2J) has to be replaced in this case with: 

(2-3) K(l,m) = ft /2 f^ II a 'M(q)- 

qGCS,, l<i<j<m 

(c) In case if is a Toeplitz kernel we noticed already that 7^+;^+; = 7^ for i < j, 
I > 1 and we denoted j n — 7o,n> Ji > 1. We conclude that di+ij+i — a^j for i < j, 
I > 1 and formula (|2.2(1 reduces to 

(2.4) K(0,n)= Y, II a ^'(P)- 

per>„ o<i<j<n 

We can compare this result with a classical formula of Verblunsky, according to 
which there exists a polynomial — V^'^i, . . . J n -U di, ■ ■ ■ , d n -i) with integer 
coeffcients so that 

n-i 

K{G,n)= ln Y\dl + V^ 

k=l 

5 



(see [B], in particular, pg. 60-61, for a comprehensive discussion of this formula). 
We see that the term j n Yik=i ^1 corresponds to the path po made of n consecutive 
rise steps followed by n consecutive fall steps. Consequently, we deduce from (|2.4(l 
that 

vw= n a ^(p)' 

p£X>„-{po}0<i<j<n 

an explicit formula that explains some of the features of V {n l □ 

3. Near tridiagonal models 

In this section we show that positive definite kernels do not have tridiagonal 
models. Instead we introduce a near tridiagonal model and then we show how 
this model is related to the set of Lukasiewicz paths. Again, in order to simplify 
the notation we consider only strictly positive definite kernels if and we assume 
if (0,0) = 1. We denote by V C Z 2 (N ) the vector space generated by the standard 
basis of / 2 (No) and we call tridiagonal model of if a family { J„} rl >o of tridiagonal 
operators (not necessarely bounded), 

&o(ra) ci(n) 
oi(n) b 1 (n) c 2 (n) 



J n = 







such that Jo 
(3-1) 



I and 



K(i,j) = e Ji 



a 2 {n) 


,J*.L. 



b 2 (n) 



J\e* , i,j>0. 



Also, in analogy with the Hankel case, we ask ak(n) > 0, k, n > 1. Since each J n is 
tridiagonal, J n V C T> so Ij3.1|l makes sense. However we have the following result. 

Theorem 3.1. There are strictly positive definite kernels with no tridiagonal model. 

Proof. We consider a strictly positive definite kernel K with 

K(0, 1) = if (0,2) = if (1,2) = 0, if (0,3)^0 

(it is easy to construct such a kernel using, for instance, Theorem l2.2ll . Let { J„}„>i 
be a tridiagonal model of if, then we deduce that 

6o(l) = d(2) = 6i(2)=0, 

which implies 
if (0,3) - e J 3 J 2 J ie * 

= &o(3) (6o(2)&o(l) + ci(2) ai (l)) + ci(3) (ai(2)6o(l) 
a contradiction showing that if has no tridiagonal model. 



6 1 (2)a 1 (l)) = 0, 
□ 



We are now trying to find a model as close as possible of being tridiagonal, which 
should reduce to (|1.4(l in case the kernel if is Hankel. Thus, we consider operators 
(not necessarely bounded) with matrix still of Hessenberg form 

&o(l) c ,i(n) co j2 (n) 



(3.2) 



Jn 



ax(l) 




&i(2) 

02(2) 





ci,2(n) 
&2(3) 



n > 1, 



with respect to the standard basis of / 2 (Nq), and with the additional conditions: 



(3.3) 



Ofc(fc) > 0, k > 1, 
Cij(l)=0, j>2, 0<i<j-l, 
Cjj (n) = 0, j > n, < i < j - 1, 
c k -i,k( n ) = a fc( fc ), A: > n, 

/ (n) = Q j- (n - 1), j < n - 1, < i < j 



We see that J n T) C £> for each n > 1. We call such a family {J„}„>i a near 
tridiagonal model of the kernel if provided that 



(3.4) K(i,j) = e JZ...J*J j ...J 1 e* , i,j>0. 

Theorem 3.2. Any strictly positive definite kernel has a near tridiagonal model. 

Proof. Let K be a kernel and K n = [K (i, j)] 0<i j <n - Then K is strictly positive 
definite if and only if K n > 0, n > 1 (as already mentioned we can assume without 
loss of generality that K(0,0) = 1). Let D n — [di.j{n)] Q<i . <n be the upper triangu- 
lar Cholesky factor of K n , therefore K n — D* n D n and <ij,i(n) > 0. The uniqueness 
of the Cholesky factor implies that 

D-n l-n+1 

d n+1 . n+ i(n + 1) 

and d n . n (k + 1) = d n>n {k) for k > n, so we can drop the label n in dij(n). We now 
construct the near tridiagonal model of K. Thus we prove by induction on n that 
there exist numbers afc(fc), bk-\{k), Cj.fc_i(fc), < i < k — 1, k <n, such that (|3.3|) 
holds and 

1 ^2,1 ^3,2^2,1 ■ ■ • Jn+l,n ■ ■ ■ ^2,1 
On n _l 0„_2 ■ ■ ■ Oq 



D 



n+1 



(3.5) 
where 



D,. 



&o(l) 
ai(l) 



Jk 



+\,h — 



6o(l) c ,i(2) 
ai(l) 61 (2) 



c ,fe_i(A;) 











02(2) 



a k (k) 



and Ofc denotes the column with k zero entries. 
For 71 = 1 we define 

6o(l) = <k>,i an d a i(l) — ^1,1 > 0; 



so that D\ 



J2,l 







Assume the statement is true up to n. We determine 



the numbers a„+i(n + 1) > 0, b n (n + 1), and Ck, n {n +1), k = 0, . . . n, n — 1, such 
that 



+ l,n+l 



— Jn+2,n+lJn+l,n ■ ■ ■ ^2,1- 



By the induction hypothesis 



Jn+l,n ■ ■ ■ JlS — 
7 



I'll 



so that we must have 



(3.6) 



ln+1 
d n +l,n+l 



X 
Xl 



cx, n {n ■ 



x n + b n {n + l)d„,„ 
a n+ i(n + l)d 

n,n 

where xq, ■ ■ ., x n are numbers uniquely determined by a,k(k) > 0, bk-i(k), and 
Ci,i(k), i < I < k < n. Since d n ^ n > 0, l|3.6|) uniquely determine the numbers 
co,n(n + 1), • ■ ■) c n _i )n (n + 1), b n (n + 1) and a n+ i(n + 1) = > o such that 

(EHfl) holds for n + 1. 

Now we can use all the numbers afe(fc), bk-i(k), Cj ; ;(fe) in order to define J n by 
(|3.2() and then (|3.5|l shows that { J n } n >i is a near tridiagonal model of K. □ 

We notice that the label n of the numbers a n (n), 6„_i(n), (^(n) is superfluous 
due to the conditions (|3.3|l . We used it in order to have a uniform definition of J„ 
in l|3.2|l but we will drop it from now on. The proof of Theorem 13.21 gives a one-to- 
one correspondence between the set of strictly positive definite kernels on No with 
K(0,0) = 1 and the set J of families of numbers bk-i, Oi,i \ k > 1, < i < I}. 
We will call these numbers the Jacobi parameters of K . In addition, we can easily 
characterize the strictly positive definite Hankcl kernels by the additional conditions 
on the Jacobi parameters: 



(3.7) 



Cij =0, / > i + 1, 
c/c,fc+i = a k +i, k>0. 

In this case the near tridiagonal model reduces to 4f> . 

The next task is to establish an explicit formula for the Cholesky factors D n in 
terms of the Jacobi parameters. First, we obtain a recursive relation for D n . 



Lemma 3.3. For n > 1, 



where 



D n = F n (l © D n _i), 



F„ 



bo 
ai 






f 








Proof. We have Dq = 1 and then 



C0,1 

bi 




C0,2 
Cl.2 



«3 



CO,: 



a i 





Do 



Assume the statement is true up to n. Then 

D 

i 





Dn 



Jn+2,n+l ■ ■ ■ ^2,1 



By the induction hypothesis, 





' D n ' 




" F n (1 © D n _i) 












and we notice that 








J,, 


f2,n- 


-1 ■ • • ^2,1 = -fn+1 



F, 



n+l 



so that 



n+l 



F 



n+l 



1 

D„_i 





J, 



71+1,71 



1 o 
£»„-i 




J2.1 



F B+ i (!©£>„). 



□ 



The matrices F n have a very simple recursive multiplicative structure. Actually 
it is convenient to make the dependence of F n on the Jacobi parameters more 
explicit and introduce 

1 bk Ck : k+l ■ ■ ■ Ck, m -1 
cik+i bk+i Ck+l,m-l 
a k+2 

bm-l 

a m 

for m > 1 and < k < m — 1. In particular, F„.o = F n . We show that the building 
blocks of F m< k (consequently, of D n ) are the 2x2 matrices 

1 b k 

a k+ i 

and the (I — k + 2) x (/ — k + 2) matrices 

1 ... c ktl 
1 



A; > 0, 



Cfe,; — 



< k < I. 



1 

Lemma 3.4. For m > 1 and < fc < m — 1, 

F mi fc = (1 © F OT! fe_|_i) G m . k , 

where 

G m ,k = Cfe,m-1 (Cfe,m-2 © 1) ■ ■ ■ (C/^fc+i © l m _fe_2) (Bfe © l TO _fc_i) . 

Proof. The proof is a straightforward calculation and can be omitted. 



□ 



Once again it is convenient to visualize all those matrix multiplications by means 
of a transmission line picture similar to the one in Figure 3. Thus, Figure 5 illus- 
trates how to calculate D3 by using Lemma 13.31 and Lemma 13.41 In particular, if 
1 is the imput at A then at B we read the expression of K(0,3) in terms of the 
Jacobi parameters, 

K(0, 3) = 6q + koaicoa + aic ,i & o + ai^ic ,i + aia 2 c , 2 . 

9 



Figure 5 suggests a connection with weighted Lukasiewicz paths. A Lukasiewicz 
path of length n is a path in the positive quadrant of the lattice Z 2 which starts at 
(0,0), ends at (n, 0), and consists of rise unit steps, horizontal unit steps, and fall 
steps of arbitrary depth. Let £„ j0 denote the set of Lukasiewicz paths of length n. 




A °i 

Figure 5. Transmission line representation for D 3 

We also consider C n ^ the set of paths of length n in the positive quadrant, 
starting at (0, 0) and consisting of the same type of steps as above, but ending at 
(n,k). We introduce a weigth on the elements of C n ^ as follows. Let p G H n ,k 
consists of steps pi, . . ., p„. Then 



fe=i 



and 



if pk is a rise step (j, I) — > (j + 1,1 + 1) for some j > 0; 
^(Pfe) = { h if Pfc is a horizontal step (j, I) — ► (j + 1, /) for some j > 0; 

if pk is a fall step (j, I) — > (j + 1, fc) for some j > Z. 




&n-l 





FIGURE 6. Passing from paths of length n — 1 to paths of lentgh n 
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Theorem 3.5. The Cholesky factor D n = [dij] 0<j - <n is given by the formula 

(3.8) d itj = W (P)' *^ (i,i)^(0,0). 

Proof. We can prove the statement by induction on n. For n < 3, (|3.8(l is seen from 
Figure 5. The general induction step is provided by Lemma l3.3l Thus, 

do,n = &0^0,n-l + Co,l<il,n-l + ■ • • + C0/n-ld n -l,n-l 
dl, n = 0<ldo,n-l + bldi^ n -i + . . . + Ci in -id n -i, n -i 

dn,n — Q>nd>n— l,n — 1) 

and these relations are precisely those obtained by passing from weighted paths of 
length n — 1 to weighted paths of length n, as showed in Figure 6. □ 

Remarks (a) For a Hankel kernel K, Theorem 13.51 reduces to well-known results 
in the combinatorial theory for orthogonal polynomials (on the real line), 0], [Hj- 
Indeed, by (|3.7|l there are no fall steps of depth other than one. In this case, the 
summation in (|3.8|l is only over Motzkin paths, which is the classical formula in 0], 
9 . It might be interesting to note that the correpsonding formula for orthogonal 
polynomials on the unit circle, (|2.4|) . involves summation over labelled configura- 
tions in Z 2 rather than over weighted paths. 

(b) There are other significant differences between the two parametrizations dis- 
cussed in this paper. For instance, the parameters {^k,j} have the following inher- 
itance property: the parameters of the kernel — [K(l,m)] 1<l are precisely 
{lk,j}i<k<j- The Jacobi parameters do not have such a property. Another dif- 
ference involves computations of determinants. Thus, we have already notice the 
formula (we assume if (0,0) = 1): 

n 

det[K(l,m)]l m=0 = l[f k J] d lv 

k=l 0<i<j<n 

and from Lemma f3. 31 we deduce 

which does not involve all the Jacobi parameters (up tu n). So the first determinant 
formula is much tighter in its parameters. 

(c) Theorem 13.51 gives . in particular, that 

(3.9) K{0,n)= W (P)> n ^ L 

pe£„, 

Since the Jacobi parameters do not have the inheritance property mentioned above, 
we cannot have formulae of this type for any K(l, n). Instead we have the following 
construction. 

For n fixed we consider admissible steps in the psitive quadrant of the following 
types: between the vertical lines x = and x — n only of Lukasiewicz type are 
allowed and between the vertical lines x = n and x = In only reflections with 
respect to the line x ~ n oi Lukasiewicz steps are allowed (and they are weighted 

li 



with the complex conjugate of the weight of the reflected Lukasiewicz step, see 
Figure 7). 




Figure 7. Admissible steps for n = 2 

Denote by K, n ^ the set of paths in the positive quadrant of Z 2 made of admissible 
steps, starting at (0, 0) and ending at (n + k, 0), < k < n. In particular, /C n! o = 
£ n fi- The weights arc defined correspondingly. With these elements, we deduce 
from Theorem 13 . 51 that 

(3.10) K{l,n)= W (P)> 0<l<n. 

(d) We notice that the proof of Theorem 12 . 21 gives a formula for the Cholesky factor 
D n in terms of Dyck type paths analogous to formula l|3.8fl ■ Q 
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